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An investigation of the longitudinal handling qualities 
of the variable stability flight simulator, converted from 
a C-llB Instrument Flight Trainer, found the qualities to 
be an unrealistic representation of aircraft motion. Non- 
linearities found in the dc servo drive circuits were 
caused by the dc servo drive motor's starting voltage, 
striction, misalinement of the motor and reduction gear 
axles, potentiometer resolution, and gearing backlash. 

The equations of motion for the T-33 aircraft were simu- 
lated with the analog portion of a hybrid computer. 
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I . INTRODUCTION 

This investigation is the third in a series of reports 
on the modification of the Link Aviation model C-llB "Jet- 
Propelled Aircraft Instrument Flying Trainer." The modif- 
ication was instituted in the Aeronautical Engineering 
Department of the Naval Postgraduate School in 1967. 

The C-llB trainer was acquired for the purpose of mod- 
ifying the trainer to a flight simulator for use as an aid 
for class demonstrations on flight characteristics and a 
research vehicle for thesis projects. The trainer was 
originally built to provide instrument training to pilots. 
It did not simulate any particular aircraft nor did it have 
the ability to change characteristics. To be useful to 
the Aeronautical Engineering Department, the simulator had 
to be capable of having its characteristics varied so 
that various aircraft could be simulated and their 
characteristics could be compared. If a simulator is 
capable of these changes it is called a variable-stability 
simulator . 

A flight simulator is a "good" simulator if it realis- 
tically represents the flight conditions that a pilot 
encounters in an aircraft. The addition of many instru- 
ments, a cockpit, aircraft noise, motion cues, and visual 
cues all aid in giving the pilot a believable flight 
experience. The additions also contribute significantly to 
the cost of the simulator. The C-llB trainer had the 
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instruments familiar to the pilot, the cockpit, and pro- 
visions for noise generation, but not the visual and motion 
cues. It was decided that the addition of visual and motion 
cues were not needed at this time for the C-llB trainer to 
become a variable-stability, fixed-base simulator. 

Lieutenant J. A. Johnson reported on his investigation 
of the feasibility of modifying the trainer in Ref. 1. In 
Ref. 2, Lieutenant Commander C. J. Sweeney described the 
conversion of the trainer to a flight simulator be replacing 
the electronic components with solid state devices and 
solving the equations of motion with a hybrid computer. The 
computer is composed of a high speed digital computer model 
9300 manufactured by Scientific Data Systems of Santa Monica, 
California and an electronic analog computer model Cl 5000 
manufactured by COMCOR, Ind. , of Anaheim, California. 

This report deals with the investigation of the longi- 
tudinal flying qualities of the simulator which are 
dependent on the electronic, mechanical, and servo elements 
of the simulator. The flying qualities of the simulator 
allow the pilot to be the closing link in the control 
loop. The pilot actuates a control in the simulator and 
the control data is changed to an electrical input for the 
computer. The computer solves the aerodynamic equations 
of motion using the control data input and the aircraft's 
characteristics and sends the solutions back to the sim- 
ulator for the deflection of the cockpit instruments. The 
pilot observes the instruments and makes the appropriate 
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control movements to perform the maneuver he desires. If 
the computer does not solve the equations correctly, or 
the instrument dials do not indicate the computer's solu- 
tions, the pilot must make additional control movements 
that are not normally required. The requirement for extra 
movements is caused by the poor flying qualities of the 
simulator and could result in the pilot forming incorrect 
opinions of the flight characteristics that are being 
tested. 

Using data for the T-33 aircraft, the flying qualities 
of the simulator were investigated in two parts: first, the 
comparison of the expected solutions to the equations of 
motion with the solutions using the analog portion of 
the hybrid computer, and second, the observation of the 
cockpit representation of the computer's solutions. The 
investigation was conducted at the Naval Postgraduate 
School, Monterey, California, during the period January 
1969 through June 1969. 
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II. COMPUTER 



During modification of the trainer, the vacuum tube 
amplifiers used for the solution of the equations of 
motion were removed with the intention of finding another 
method to solve the equations. Instead of designing and 
building solid state amplifiers to fit into the simulator, 
it was economically and operationally feasible to use the 
hybrid computer facilities of the Naval Postgraduate 
School to solve the motion equations. The computer also 
offered a very rapid method of varying the characteristics 
of an aircraft, either manually or by a precomputed 
program. The instrument calibration factors for convert- 
ing voltages to dial deflections were easily accomplished 
on the computer. 

A. FORCE AND MOMENT EQUATIONS 

Four axes systems are commonly used in writing the 
equations of motion: inertia, body, wind, and stability. 
Each system is a right-hand Cartesian coordinate system 
with all but the inertia system having its origin at the 
center of gravity, c.g., of the aircraft. The inertia 
system is fixed in inertia space having its origin on the 
earth's surface and the x-axis fixed in the direction of 
true north. The body axes have the x- and z- axes fixed 
in the plane of symmetry and the x-axis in the direction 
of the aircraft's nose. The wind axes have the x- and 
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z-plane rotated to have the x-axis in the direction of the 
velocity vector. The stability axes are between the body 
and wind axes with the x-stability axis the projection of 
the x-wind axis into the plane of symmetry. The rotation 
angle of the x-body axis to the x-stability axis is the 
angle of attack, a, and the rotation angle of the 
x-stability axis to the x-wind axis is the sideslip angle, 

3 . The body, stability, and wind axes are related to the 
inertia axes system with the Euler angles of pitch, 6, 
roll, 4>, and yaw, ip. Body and stability axes relationships 
are shown in Figure 1 and the axes relationship with the 
Euler angles are shown in Figure 2. 

Three of the four axes systems were utilized for the 
longitudinal equations in the computer: the body axes for 
the moment equation and the Euler angles, the stability 
axes for the force equations , and the inertia axes for 
the altitude and rate of climb. The full equations 
[Ref. 3] are: 

Moment Equation About y-Body Axis 

M, = I q. + (I - I ) p. r, - I ” Pu) 

b y^b ' X z *^b b xz b 

Force Equations in x-and z-Stability Axes 

Xg=m (r^^cos a -pj^sin a)+ g (sin 6j^ cos a-cos 6j^cos <}>|^sin a) ] 

- T cos (a + a^) 

Z = m [v (p, cos a + r, sin a) - u„ (q. - a) 

s s ^b b s ^b 

- g(cos 9j^ cos ()>j^ cos a + sin 6j^ sin a) ] + T sin (a + a^) 
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FIGURE I. RELATION OF BODY AND STABILITY 

AXES 




FIGURE 2. MASTER DIAGRAM OF EULER ANGLES 
AND AXIS SYSTEMS USED IN SIMULATION 
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Euler Angles in x - 2 Body Axes Plane 
6b = qb cos 0b - Tb sin c()b 



Height in z-Inertia Axis 

-h = -u cos a sin 6. + v sin d) cos 6 + u 
s D s ^ s 



sin a cos (|) cos 



6 



In order to simplify the above equations, the following 
assumptions were made; 

1) Motion in x-z plane only, <}> = \p = ^ v = 0 

s 

2) Assume a and 6 are small, so that 

sin a s a cos a = cos 0 s cos (a + a^) = 1 
sin 0=0 

3) Products of a, 0, p, and r are negligible. 

4) Assume u = V and V = V + /Vdt 

s o 

5) Thrust component T sin (a + a^) is negligible. 

Use of these assumptions, which were made in order to use 
small perturbation theory and three degrees of freedom, 
resulted in the following simplified equations that were 
used in the simulation: 

Mu = I 4u 

b y ^b 

Xg = m [V + g (0b-a) ] - T 
= m [-V^ -gl 

' ‘’b 

-h = -V (0b-a) 

The first derivative of the varying quantities can be 
obtained by rearranging the simplified equations; 
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a 





h = V (6, - a) 

b 

The above equations are the general linearized equations 
for any aircraft. To simulate a specific aircraft, expres- 
sions with the flight characteristics of the aircraft are 
substituted for X , Z , and M, . 

S S iD 

B. AERODYNAMIC EQUATIONS 

Because the Lockheed T-33 aircraft is used for Combat 
Readiness Training Flights at the Naval Postgraduate 
School, it was decided to simulate the T-33 on the modified 
trainer using data from Ref. 3. The aerodynamic force and 
moment equations are: 



+Cp a+Cp 

o a 



a^)-q^S^(Cp ^V+Cj^ aAV+c' a^AV) 



O o a 
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-■Jo^ <^L ■"'^L ‘^e' “*''> 

o a 6 o o a 

£ 

E 

I (M a+M- a +M q + M~ 6„) 
y a a q ^ ^ 

Analog diagrams of the equations for simulation are shown 
in Figures 3 , 4 , and 5 with the potentiometer settings 
listed in Table I. 

The longitudinal mode was excited with a step input of 
one degree elevator deflection to test the data, equations, 
and the analog setup. Results of the positive and negative 
step functions are shown in Figures 6 and 7. The unreal- 
istic rate of climb indications of 6000 feet per minute 
after five seconds of elevator deflection and 11,000 to 
13,000 feet per minute after ten seconds of deflection 
suggested an error in the program that could not be found. 
The phugoid period of 65 seconds observed in the 6 response 
for the negative step input was an acceptable period. The 
phugoid period for the positive input could not be observed 
because the computer's amplifiers overloaded after 30 
seconds of input. Due to the lack of the actual known 
results, the obtained results could not be compared to 
determine the cause of the error. A correct analog setup 
was verified by comparing the hybrid computer results 
with results from an IBM 360 computer. Using the same data 
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FIGURE 3. ANALOG SIMULATION CIRCUIT (PART I) 
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FIGURE 4. ANALOG SIMULATION CIRCUIT (PART 2) 
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FIGURE 5. ANALOG SIMULATION CIRCUIT (PART 3) 
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TABLE I 



ANALOG POTENTIOMETER SETTINGS 



P o ten ti ome ter 


Reference 3 


Reference 5 


D^/2000 


.5100 


.4740 


D /200 
a 


.7050 


.3635 


C„2/20 


.0408 


.2475 




.1382 


.8100 




.2389 


.1190 




.1727 


.1550 


KTS 


.5921 


.5921 


1/t 


.1000 


.1000 


1/5 


.2000 


.2000 


Jiig/sv , 300 


.2164 


.2241 


V /lOOO 
o 


.6100 


.6130 


100/m 


.2597 


.2507 


L^/100,000 


.1240 


.1284 


L /10,000 
a 


.5302 


.5650 




.1739 


.1849 


L^v/100 


.4065 


.2955 


L /lOOO 


.3016 


.3016 




.1000 


.1000 


“«E V^“°° 


.9523 


.9580 


M„Iy/5000 


.4877 


.7200 


M-I /500 
a y ' 


.5282 


.2050 


M I /500 
q y 


.5318 


.4530 


2865/Iy 


.1384 


.1384 


mg/50 ,000 


.2480 


.2568 


57.3/mV^ 


.1221 


.1172 


vyS7.3g 


. 3310 


.3320 


1/100 


.0100 


.0100 


10/57.3 


.1745 


.1745 


1/2.5 


.4000 


.4000 
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FIGURE 6 . RESPONSE TO A POSITIVE ELEVATOR INPUT (Ref. 3 Data) 
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and equations in the two computers, the similar results 
obtained from both of the computers indicated no error 
in the analog setup. 

Data from Refs. 3, 5, and 6 was compared to determine 
if possible errors in the input data to the computer would 
cause the obtained results. Minor variations in the 
initial conditions between the references were not con- 
sidered significant enough to alter the results. Although 
most of the data in Refs. 3 and 5 was different, the results 
shown in Figures 8 and 9, using Ref. 5 data, are not a sig- 
nificant change from the results with data from Ref. 3 
shown in Figures 6 and 7. The data from Ref. 5 is listed 
in Table I. The rate of climb indications were 3,600 feet 
per minute after five seconds of elevator deflection and 
9,000 feet per minute after ten seconds of input. The 
phugoid period increased to 78 seconds for the negative 
step input, which is an acceptable response. 

In order to use the computer's solutions to the equa- 
tions of motion for checking the cockpit's instrument move- 
ment, the aircraft's characteristics were modified until 
the positive and negative inputs did not overload the 
computer. In some cases, the computer solution signal to 
an instrument was small, so it had to be replaced with a 
large amplitude sinusoidal signal in order to obtain large 
instrument dial deflections. 
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III. COCKPIT 



The C-llB trainer was modified in two ways: the physical 
dimensions of the trainer were shortened, and the vacuum 
tube amplifiers and motor-generators were replaced with 
solid state devices and dc motors, respectively. The 
interior of the cockpit remained basically unchanged. The 
trainer's instruments were retained because of their simi- 
larity to present aircraft instruments and the pilot's 
familiarity with them. The navigation equipment was retained 
for cockpit realism but was not made operational. 

A. CONTROLS 

The cockpit of the simulator represents a typical 
cockpit and has the normal control devices found in most 
cockpits: control stick for aileron and elevator deflec- 
tion, foot pedals for rudder deflection, throttle, speed 
brakes, landing gear, and flaps. The speed brakes, landing 
gear, and flaps were not connected but provisions were 
made for their eventual operation. The shortening of the 
trainer required the linkage for the stick and rudder 
controls to be altered and relocated in the simulator. 
Because the elevator portion of the control stick and 
throttle are the only controls that affect the simulator 
in the longitudinal mode, they were the only ones considered 
in this report. 
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The movement of the linkage for the elevator control 
had no adverse effect on the force required to control the 
simulator. The stick force that the pilot must overcome 
is regulated by a dynamic pressure servo. The servo is 
presently manually set at a predetermined value but can 
be programed to vary with altitude and true airspeed. 

All control force variation due to linkage alteration was 
absorbed in the dynamic pressure servo setting. 

The amount of elevator deflection used in the equations 
of motion is determined from the stick deflection and the 
elevator deflection potentiometer at the computer. The 
deflection of the control stick moves the arm of a tapped 
potentiometer on the elevator control mechanism and the 
electric signal is sent to the elevator deflection poten- 
tiometer. This potentiometer controls the elevator 
effectiveness of the simulator by allowing the amount of 
elevator deflection per unit of stick deflection to be 
changed. A plot of the signal from the tapped poten- 
tiometer versus the stick deflection is shown in Figure 10. 
The nonlinear portion of the figure was caused by the 
mechanical linkage in the elevator control mechanism. 

Because of the improper instrument dial response to a 
signal from the computer, that will be discussed later, 
the effect of the stick movement in the equations of 
motion could not be determined. 

The throttle controls the amount of thrust delivered 
by the simulated engine. Because the engine characteristics 
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FIGURE 10. ELEVATOR CONTROL STICK CHARACTERISTICS 



were not known, the engine was simulated on the computer 

as having a thrust of 1020 pounds at an altitude of 

23,000 feet [Ref. 3]. The engine simulation included a 

first-order approximation with a one second time constant. 

The expression for the engine simulation was 

T = (1-e T , • 4 . (where T is the throttle 

computer cockpit 

signal and t is the time constant) . The effect of throttle 
variation was not investigated. A more complete simula- 
tion could be accomplished using the digital part of the 
hybrid when sufficient data is available. 



B. INSTRUMENTS 

The instruments used in the simulator are typical 
instruments found in present day aircraft. The driving 
circuits of the instruments were modified to operate from 
the dc computer signal by replacing the ac amplifiers and 
motor-generators with dc motors and solid state amplifiers . 
Modification of the following instruments was described 
in Ref. 2: 



accelerometer 
airspeed (indicated) 
altimeter 
angle of attack 
gyro horizon (bank) 
heading (RMI) 
needle and ball 
tachometer 
vertical speed 

The gyro horizon (pitch) could not be modified to 
operate from the dc computer signal, so the dc pitch signal 
was modified to produce an ac signal. The general 
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operation of the modified circuit is to attenuate the dc 
computer signal prior to entrance into the amplifier sec- 
tion. The amplified signal passes through either a 
P-N-P or a N-P-N transistor and drives a dc motor in the 
proper direction. The motor, through gear trains, moves 
the follow-up potentiometer and the synchro transmitter 
for the instrument. The motor is driven until sufficient 
voltage from the follow-up potentiometer is directed back 
to the amplifier to reduce the amplified signal to zero. 

The synchro receiver in the instrument receives the signal 
from the synchro transmitter and moves the dial of the 
instrument. A typical modified circuit is represented 
by the Modified Airspeed Indicator Drive Circuit shown in 
Figure 11 [Ref. 2] . 

The accelerometer , airspeed (indicated), altimeter, 
angle of attack, vertical speed, and gyro horizon (pitch) 
instruments were utilized during the longitudinal mode 
investigation. Because of the lack of valid engine infor- 
mation, the tachometer was not considered. The initial 
conditions for the instruments and their calibration are 
shown in Table II. 

The motion of the instruments , when a signal was 
applied, was observed to be erratic and jerky. The non- 
linear motion was determined to be caused by the dc motor 
and the motion of the arm of the follow-up potentiometer. 
Linear relationships existed between the input signal/ 
attenuated signal and between the synchro transmitter 
signal/instrument dial. Because the synchro transmitter is 
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FIGURE II. MODIFIED AIRSPEED INDICATOR DRIVE CIRCUIT 



TABLE II 



INITIAL CONDITIONS AND CALIBRATION 
OF FLIGHT INSTRUMENTS 



Instrument 


Initial 

Condition 


Volts 

D.C. 


Calibration 




Accelerometer 


3.2 g's 


+ 100 
- 100 


-0.8 g's 
+6.5 g's 




Airspeed 


250 kts. 


+ 100 
- 100 


92 kts. 
405 kts. 




Altimeter 


23,000 ft. 


+ 100 
- 100 


20,030 ft. 
26,020 ft. 




Angle of Attack 


12 deg. 


+ 100 
- 100 


26 deg. 
-1 deg. 




Gyro Horizon 
(pitch) 


2 deg. up 


+ 100 
- 100 


65 deg. up 
45 deg. down 




Vertical Speed 


0 


+ 100 
- 100 


4000 ft. /min. 
4200 ft. /min. 


down 

up 
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fc© th© follow-up potentiometer arm, 

V§itSf§ th@ f©il©W“Up potentiometer arm can 

the iial defleetion of the instrument. An 
§M§fiipie §1 a fe^pigal eireuit'g nonlinear characteristics 
i§ §h©wfl the airepeed eireuit characteristics in 
Figute ii. 

¥he de m©t©r did n©t react immediately to a signal 
fr©ffl the amplifier heeause ©f the motor's excitation 
voltage f §ti©ti©n, and misaiinement of the motor and 
red\i©t4©n ^ear axlea t a further delay was noticed in the 
inatrument dial due t© gearing backlash. A smaller voltage 
was required t© start the m©t©r rotating when the motor was 
h©t attaehed t© the redueti©n gears than when the motor 
was attaehedi S'he voltage required to start the motor 
ai§© varied from instruiftent to instrument because of the 
degree ©f miaaiinement ©f the motor and reduction gear 
axles V the nuifeer ©f gears in the gear trains between the 
®(©%©r and the serv© transmitter and between the motor and 
the f©ll©w=-np potentiometer-, and the sise of the motor used, 
two sioes of d© motors were used because of insufficient 
amounts ©t one sises 

the motor ''s nonlinearity was affected by the resolution 
of the wire -wound follo’^w=-up> potentiometer [[®ef , , ihe 

voltage range of bhe ^ollow^-up control was p>lus IS "^its 
to minus 1§ volts-, h. siftall movement of the arm ®)r©du! 2 ed 
a large feediaeh signal in oonparison to the irspat signal.. 
Reversal ot the motor was delayed luntil sufficient voltage 
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FIGURE 12. AIRSPEED CIRCUIT CHARACTERISTICS 



to overcome the backlash was directed to the amplifier by 
either the input or potentiometer signals. Figure 13 
shows the problem of resolution and backlash for a sinus- 
oidal input. A tachometer feedback from the motor was 
not feasible because the motor's rotation was small once 
the backlash of the circuit was exceeded. 

A reduction in the nonlinearity was attempted by in- 
creasing the input to the amplifier, increasing the ampli- 
fier gain, adding a bias to the signal, and increasing 
the follow-up gain. Increasing the ratio of the attenuat- 
ing resistor to the total input signal resistance resulted 
in increasing the signal to the amplifier for the same 
computer input signal. Although this ratio change of 
resistors resulted in larger movements of the motor shaft, 
the change was not used because it also caused larger 
instrument dial deflection. Reducing the amount of instru- 
ment dial movement in order to observe incremental dial 
deflections resulted in the size of the computer input 
signal being reduced. If the signal was reduced too much, 
the probability of signal error due to noise was increased. 

Increasing the amplifier gain by changing the feedback 
resistor resulted in the same problems encountered with the 
increased input signal plus an additional problem by the 
amplification of the resolution error in the follow-up 
potentiometer. The addition of a bias voltage was bene- 
ficial when movement was restricted to a single direction 
but hindered the movement of the motor during motor 
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FIGURE 13. COMPARISON OF AIRSPEED SERVO INPUT AND FOLLOW-UP POTENTIOMETER SIGNALS 



reversal. An increase in the follow-up potentiometer gain 
amplified the resolution error and increased the staircase 
effect of the potentiometer. 

The following corrective actions were considered but 
not implemented because of the lack of time: circuit 
redesign, decrease of the voltage range of the follow-up 
potentiometer, and computation of nonlinear describing 
functions to use in the simulation. The circuit redesign 
would present a problem in matching the forward and reverse 
driving voltages to the motor. With the redesign of the 
gear ratios between the motor/synchro transmitter and the 
motor/follow-up potentiometer arm, a larger motor shaft 
rotation could produce a small movement of the instrument 
dial's synchro transmitter and the follow-up potentiometer 
arm. This would allow a larger signal to the amplifier for 
a desired dial deflection. 

Decreasing the voltage range on the follow-up poten- 
tiometer would require a larger movement of the follow-up 
arm in order to feedback a sufficient voltage to cancel 
the input signal. The larger movement would also cause 
the instrument dial to have a large movement because of 
the gearing. To obtain the small movements required in 
the dial deflection, smaller input signals would have to 
be used which have the problems of signal noise error. 

The describing function method, in order to take into 
account the motor's starting voltage, stiction and the 
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backlash, would be very involved because each circuit has 
different characteristics. 

The nonlinearities might be reduced by acquiring 
different dc motors that do not require a large starting 
voltage . 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



The longitudinal flying qualities of the modified C-llB 
Instrument Flight Trainer are unsatisfactory and do not 
represent the flying qualities of a present day aircraft. 
The nonlinear instrument dial movement interferes with 
the pilot's control of the simulator by requiring extra 
control movements for a desired maneuver. The lack of 
accurate instrument indications results in erroneous con- 
clusions by the pilot of the flight characteristics being 
simulated. Because the lateral direction instruments 
utilze similar servo circuits as the longitudinal mode, 
the lateral mode can be assumed to be unsatisfactory, too. 

The unrealistic flight data results for the T-33 
aircraft requires a "hit or miss" guess at the correct data 
to be used in the simulation. The conversion of the 
flight trainer to receive a dc signal from the hybrid 
computer resulted in a simulator whose aircraft stability 
derivatives are easily changed. 

The versatility of a variable stability simulator for 
use in demonstrating flight characteristics to academic 
classes and for providing a research vehicle for the in- 
vestigation of flying qualities of aircraft, pilot response 
techniques, new flight instrumentation, and variable- 
control systems justify the continuation of the modifi- 
cation of the flight trainer. Improvements in the 
instrument dial representation could be accomplished with 
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acquisition of low torque dc servo motors, do potentiometer 
instruments, or redesign of the present servo motor 
amplifier circuit. Reliable aircraft data for the T-33 or 
any single engine jet aircraft is required for proper simu- 
lation. The expansion of the equations to six-degrees of 
freedom with their nonlinear terms is possible with the 
hybrid computer. 

Preliminary indications are that the instabilities 

of the computer solution lie in the nonlinear terms, 

2 

such as a and aAV. It is recommended that further work 
in this area be undertaken in two steps; first, eliminate 
the nonlinear terms and examine the results of the com- 
puter solution, and second, investigate the reasons for 
the unsatisfactory results when the nonlinear terms are 
used. 
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